1986 


N87- 


^ 37 - 3 

16779 3li 




NASA/AseEj SUMMER FACULTY 


r»r»r r AMOU T 
i: uuuvnuui. 


D 


DDA^DaM 


MARSHALL SPACE FLIGHT CENTER 
THE UNIVERSITY OF ALABAMA 


SHARPLY CURVED TURN AROUND DUCT FLOW PREDICTIONS 
USING SPECTRAL PARTITIONING OF THE TURBULENT 
KINETIC ENERGY AND A PRESSURE MODIFIED WALL LAW 


Prepared By: 

Academic Rank: 

University and Department: 

NASA/MSFC : 

Laboratory: 

Division: 

Branch: 

MSFC Counterpart: 

Date: 

Contract No: 


Assistant Professor 

Memphis State University 
Mechanical Engineering 


Systems Dynamics 
Atmospheric Sciences 
Computational Fluid Dynamics 

Charles Schafer 

August 28, 1986 

NGT 01-002-099 

The University of Alabama 


XXXVII 


SHARPLY CURVED TURN AROUND DUCT FLOW PREDICTIONS 
USING SPECTRAL PARTITIONING OF THE TURBULENT 
KINETIC ENERGY AND A PRESSURE MODIFIED WALL LAW 


BY 


L. Michael Santi 

Assistant Professor of Mechanical Engineering 
Memphis State University 
Memphis, Tennessee 


ABSTRACT 


Computational predictions of turbulent flow in sharply curved 180 
degree turn around ducts are presented. The CN52D computer code 
developed at Marshall Space Flight Center is used to solve the 
equations of motion for two-dimensional incompressible flows 
transformed to a nonorthogona '» body-fitted coordinate system. This 
procedure incorporates the pressure velocity correction algorithm 
SIMPLE— C to iteratively solve a discretized form of the transformed 
equations. A multiple scale turbulence model based on simplified 
spectral partitioning is employed to obtain closure. Flow field 
predictions utilizing the multiple scale model are compared to 
features predicted by the traditional single scale k-e model. Tuning 
parameter sensitivities of the multiple scale model applied to turn 
around duct flows are also determined. 

In addition, a wall function approach based on a wall law 
suitable for incompressible turbulent boundary layers under strong 
adverse pressure gradients is tested. Turn around duct flow 
characteristics utilizing this modified wall law are presented and 
compared to results based on a standard wall treatment. 
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NOMENCLATURE 


DIFINITION 


production constant in ss k-e turbulence model 
(=1.44, see equation (5) ) 

dissipation rate constant in ss k-e turbulence model 
(=1.92, see equation (5) ) 

production constant in e p transport relation 
(=1.6, see equation (9) j 

dissipation rate constant in e p transport equation 
(see equation (9) ) 

spectral transport constant in e t transport equation 
(=1.15, see equation (10) ) 

dissipation rate constant in e t transport equation 
(see equation (10) ) 

diffusivity constant in two-equation models 
(=0.09, see equation (3) ) 
duct width 

turbulent kinetic energy 

large eddy turbulent kinetic energy 
(see Figure 2) 

intermediate spectral range turbulent kinetic energy 
(see Figure 2) 

large eddy upper wave number limit 
(see Figure 2) 

dissipative range lower wave number limit 
(see Figure 2) 
inlet plane pressure 

turbulent energy production term 
(see equation (6) ) 

TAD centerline radius of curvature 

time 

x-directed velocity component, or 
TAD longitudinal velocity component 
dimensionless universal velocity 

Reynolds stress component 
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DEFINITION 


time average longitudinal velocity 

duct average longitudinal velocity 

y directed velocity component, or 
TAD spanwise velocity component 
friction velocity 

Cartesian coordinate direction 

Cartesian coordinate direction 

dimensionless y (=yv*/V ) 

turbulent kinetic energy dissipation rate 

large eddy dissipation in the cascade 

intermediate range eddy dissipation in the cascade 

spanwise coordinate in BFC transformed system 

angle measured from TAD inlet plane 

absolute viscosity 

effective total viscosity (=y+u t ) 

effective turbulent viscosity in Boussinesq 

approximation 

kinematic viscosity 

kinematic turbulent viscosity 

longitudinal coordinate in BFC transformed system 
density 

large eddy turbulent kinetic energy Prandtl number 

intermediate range eddy kinetic energy Prandtl 
number 

large eddy dissipation rate Prandtl number 
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SYMBOL 

0e t 


DEFINITION 

intermediate range eddy dissipation rate Prandtl 
number 

stream function 


Re 
CFD 
BFC 
SIP 
SSME 
TDMA 
ms k- £ 
ss k-e 

mwf 

* 


Reynolds number 

computational fluid dynamics 

body-fitted coordinates 

Stone's implicit procedure 

space shuttle main engine 

tridiagonal matrix solver 

multiple-scale k-e model 

conventional single-scale k-e model 

pressure modified wall function 

superscript denoting dimensionless quantity 
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I 

INTRODUCTION 


The limitations of current "state of the art" turbulence models 
are quickly realized when applied in conjunction with computational 
fluid dynamic (CFD) analysis of typical turbomachine flow 
configurations. The tortuous nature of the flow path as well as the 
occurrence of complex phase and composition changes present obstacles 
that preclude totally satisfactory modeling. 

Much attention has been directed toward turbulence model 
improvements based upon the incorporation of additional pressure 
strain effects and/or wall contributions in streamfields with 
significant curvature [1-9]. Due to the computational complexity of 
CFD analyses applied to typical turbomachine configurations, and to 
the inherent limitations of the popular k-e model upon which most 
curvature modified constitutive relations are based, there is no 
general consensus on the desirability of any specific curvature 
modified model. 

Several numerical methods utilizing body-fitted coordinates have 
been developed which are suitable for solving the incompressible 
Navier-Stokes equations in complex passageways [10-14]. Each method 
employs a specific strategy for determining the pressure field such 
that the flow field is divergence free. Although there are many 
heuristics associated with these CFD techniques, it is apparent that 
sufficiently robust computational strategies are now available for 
analysis of complex, curvature-dominated flow regimes. It is, 
therefore, obvious that accuracy of turbulent flow prediction depends 
largely on the quality of the turbulence model. 

Recently, several turbulence models have been developed [15-18] 
which address the most serious drawback of two-equation closure 
models, namely, the characterization of the w.hole spectrum of 
turbulent motion by a single set of scales. This simplification 
inherently assumes equilibrium spectral energy transfer and can be 
appropriately applied only to cases where the mean flow evolves 
slowly. This type of flow condition almost never exists in typical 
turbomachine applications. The modest successes of two— equation 
models in predicting complex curvature dominated flows is more a 
testimony to the perseverance of researchers in this arena and the 
demand driven development of computational heuristics to accomodate 
experimental observation than to any inherent physical applicability. 


XXXVII-1 


The so-called "multiple-scale" turbulence models feature a two 
regime partition of the energy spectrum. The turbulence properties 
associated with large scale energetic vortices are related to mean 
velocity gradients while the turbulence dissipation occurring 
primarily in small scale eddies is related to energy transfer by its 
own action rather than to the overall kinetic energy. 

In section III df this report, a brief outline of a specific 
multiple-scale model developed in reference [18] is presented. This 
turbulence model, as implemented in the SIMPLE-C based computer code 
[14] developed at Marshall Space Flight Center, was applied to analyze 
high rate turbulent flows in the two-dimensional turn around duct 
(TAD) configuration described in Figure 1. The sharply curved section 
of this geometry is characteristic of the space shuttle main engine 
(SSME) fuel side turbopump TAD. In addition, this configuration is 
reasonably representative of general, curvature dominated, internal 
flows and as such provides a meaningful test case for measuring the 
predictive capability of the specific multiple-scale turbulence model 
employed. 

Computationally derived TAD flow predictions generated using the 
CNS2D implementation of the twin scale turbulence model are compared 
to results obtained using the standard k-e model and to results 
obtained using a curvature modified wall function [9] in conjunction 
with the traditional k-e model in section IV. In all analyses, the 
widely used wall function approach is used for the treatment of near 
wall boundary conditions [19,20] in order to achieve computational 
efficiency. 
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Figure 1. Geometric configuration for TAD analyses. 


II 

OBJECTIVES 


The general objectives of this investigation are threefold. 
Contributions to the data base of computational results involving 
curvature dominated flows are to be provided. Secondly, comparisons 
of computationally derived results utilizing several turbulence model 
modifications are to be cataloged in order to test the sensitivity of 
flow prediction to turbulence model. Finally, a tentative evaluation 
of parent CFD code reliability and robustness is sought by 
quantitative comparison with experimental results and by qualitative 
comparison with previous computational results and extrapolated 
experimentation. 


The specific objectives of this study are listed below. 

1) To provide computational predictions of TAD flow field 
characteristics at high flow rates utilizing the CNS2D computer code. 


multiple-scale 

implementation. 


C ^ 1 J 
JL. AC XU 


t ur bu 1 ence 
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model within the CNS2D 


.. M 4 M M J- U _ 
uaiii^ Luc 

computer 


3) To study the effects of a modified wall function treatment on 
TAD flow predictions at very large flow rates. 

4) To compare the phenomenological characteristics of TAD flows 
provided by the multiple-scale turbulence model with the predictions 
of traditional two-equation models. 

5) To estimate tuning parameter sensitivity of the specific 
multiple-scale turbulence model employed. 
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Ill 

BACKGROUND 


The CNS2D computer code was developed by Y. S. Chen at NASA's 
Marshall Space Flight Center. It provides a numerical procedure for 
solving the discretized equations of motion for two-dimensional flows 
utilizing nonor thogonal body-fitted coordinate (BFC) systems. The 
equations of motion are transformed to a curvilinear coordinate system 
appropriate to the geometric configuration under analysis. For 
complex geometries, individual cells defined by the spatial grid 
discretization of the problem are locally transformed into regular 0-1 
square volumes in (£,n) space by means of a bilinear transformation. 
For incompressible Newtonian fluids, the motion governing equations in 
Cartesian coordinates can be written as 


+ F y = S 


(1) 


where 


E = 



r? v 

- -^U.VV-'-Uy 


vv 


-A v y 


Co 

s -iit la. a N + r„ i 

1 7 *'x / * ■ 7 - r 

I f * i 1 1 \ . f . . * x _ O I 

y/y * yj 

In the (£,n) transform domain these equations become 


E S^x + E n n x + F C^y + F u ri y = s ^' n) • 

The metrics C x , n„, C y , and U are computed numerically using second 
order central differencing. 


For turbulent flow calculations, the molecular viscosity in 
these equations is replaced by an effective viscosity h e f f = U + Pf- / 
where ^ is an effective turbulent viscosity introduced via the 
Boussinesq approximation. The turbulent viscosity must be supplied by 
appealing to an appropriate turbulence model. In the popular 
model, u t is related to the turbulent kinetic energy, k, and its 
dissipation rate, e , by the expression [19] 

* p Cjjk 2 /s . (3) 

where k and e must themselves satisfy differential transport equations 
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of the form 


(pO^k - U e ff k x /°k } x * P ( p r “ e > 

i i 

(pUi - Ueff^Tc /°e ) x 88 P< c i p r ~ C 2 e)e/k (5) 

i i 

In the above equation the kinetic energy production term is given by 
the relation 


P r * CJj (k 2 /e ) [(u y +v x ) 2 + 2 (u x 2 +v y 2 ) ] (6) 

and the model constants are generally accepted to be approximately 

C u = 0.09 o k = 1.00 o e = 1.30 

C 1 = 1.44 C 2 = 1.92 

based on experimentation. 

Solution of the (£,r|) transformed version of equations (1),(4), 
and (5), in discretized form, is obtained iteratively by using the 
pressure velocity correction algorithm, SIMPLE-C [21]. Second order 
upwinding is used to approximate the convective terms. A grid 
staggering scheme, described in detail in reference [14], is employed 
to remove the difficulty in solving the Poisson pressure correction 
equation which couples the pressure and velocity fields. The system 
of algebraic equations arising from the discretization process is 
solved using a combination of the tridiagonal matrix solver (TDM A) 
[22] and Stone's implicit procedure (SIP) [23]. The wall boundary 
approach is used for the treatment of near wall boundary conditions. 

If the inherent computational limitations of the well tested 
SIMPLE based algorithmic approach are accepted as a necessary price 
for meaningful approximation of complex fluid flows, then improvements 
in predictive capability must come either from improved turbulence 
models or from modified wall boundary treatments. In the context of 
two-equation turbulence closure models, a single time scale is used to 
characterize all turbulent motions. Because turbulence comprises 
fluctuating motions with a spectrum of sizes and time scales, it is 
obvious that mathematical models providing a single time scale are 
fundamentally limited. Therefore, the motivation for multiple-scale 
modeling is readily apparent. 

The theoretical basis for the multiple-scale ir.cdel employed vi 
this study is described in detail by Hanjalic et al. [16]. A tyv ica^ 
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energy spectrum for shear flow at high Reynolds numbers is displayed 
in Figure 2. The quantities K L and K 2 in this figure denote 
respectively the wave number above which no significant mean strain 
production occurs and the largest wave number at which the viscous 
dissipation of turbulence energy is unimportant. Energy leaves the 

"proauct i on" region ^ rv ^ r^/ au « — — ; * - - 

region (K>K 9 ) at a rate e. Between thise regions is an intermediate 
range of wave numbers in the so-called "transfer region with a 
representative spectral energy transfer rate 


Although the description provided by Figure 2 is simplistic, it 
does provide for zonal shape changes in the energy spectrum whereas 
traditional single-scale models inherently assume spectral equilibrium 
with e s £. ■ e • Recognizing that sharply curved internal flows are 

suscef&ible to significant recirculation with a high degree of 
spectral nonequilibrium, the potential for turbulence model 
improvement using a multiple-scale concept is apparent. 


In order to utilize the multiple— scale concept in a manner 
analagous to the two-equation model approach, it is necessary to 
develop transport equations for kp, k^., e p , and e^. It is assumed 
that spectral equilibrium exists “between the transfer region and 
dissipative region, i.e. e = e^.. 


For practical applications the gradient formulation suggested by 
Launder and Spalding [19] was employed by Chen [18]. The model 
transport equations are thus given by the relations 


r\ k 



U>* 







# » 
i / ) 


( 8 ) 


De f 

Ot 




(C„p r -c^s) 


(9) 
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WAVE NUMBER K 


Figure 2. Energy Spectrum for multiple scale k-e model. 
See reference [16]. 
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Dt 





(c«,S- 



( 10 ) 


In equations (7) through (10), the first term on the right hand 
side obviously represents the flux of the respective turbulence 
property by diffusion* is the rate of production of kinetic energy 

by mean velocity gradients (see equation (6) ), e p is the energy 
transfer rate from the large scale disturbances to smaller scale 
eddies in the cascade, and e t is the flux of energy through the 
"transfer" region into the "dissipation" range. The time scale kp/e p 
used in equations (9) and (10) has the property that when e = e^ =e , 
the multiple— scale model recovers the single— scale k— e model [18j. 

Development of the model coefficients appearing in equations (7) 
through (10) is detailed in reference [18]. These are given below. 


°k 


P 



1 


P 


a~ = 1.22 

e t 


-pi 


= 1.6 


'P2 


= 1.8 - 


0.3 [l-k t /k p j/[l+k t /k p ] 


r, . 

'Cl 


= 1.15 


' Cl 


= 1.8 


c- yj 


:atio k t /k p , is especially 
between intermediate range 


The value of C V 2 , dependent upon the ratio 

important for it is here that cross-talk between intermediate rang* 
disturbances and the large scale eddies is introduced. If Cp^ were 
chosen to be a constant, the transport equation for e would be 
largely independent of the intermediate scale motions since neither e^ 
nor would explicitly appear. The form of C p 2 given above was 
chosen by Hanjalic et al. [16] and modified by crfen [18] in order to 
best agree with a variety of test flow data. In the context of the 
gradient diffusion Boussinesq type approximation, the suggested form 
of the eddy viscosity is expressed as [18] 


v t * Si (k P + k t ) k P /E P 


( 11 ) 
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where kp/ep is used in place of k/e as the appropriate time scale. 

It should be noted that a ddi tional production terms due to extra 
strain rate appear in the DuTuj/Dt term of the Reynolds stress 
transport equations when body-fitxed coordinates are employed in the 
analysis of curved shear flows [5,27]. These terms cannot be 
neglected and traditionally are included in the energy production rate 
term. In the (£,n) transformation domain, the rate of kinetic energy 
production would then be represented by the transformed version of 
equation (6) plus production terms due to extra strain which properly 
arise in the expression for Dk/Dt. This heuristic production increase 
is not included in the CNS2D computer code. The formal transformation 
represented by equation (2) is rigorously applied to all transport 
relations. The effect of the absence of these so-called extra strain 
production terms on the computational results that follow cannot be 
assessed without further investigation. 

The multiple-scale turbulence model provides a particularly 
simple method for potential improvement in the computational 
prediction of turbulent flows by incorporating a degree of spectral 
information. However, improper application of conventional wall 
functions to turbulent boundary layers under adverse pressure 
gradients contributes to the discrepancies observed between 
computational predictions and experimental measurements. Despite 
this, the logarithmic law velocity profile remains the conventional 
method of providing boundary conditions for the governing transport 
equations of turbulent motion. The inadequacy of this profile 
assumption for increasingly large adverse pressure gradients is 
graphically depicted in Figure 3. In the context of a wall functxon 
approach it is necessary to utilize a wall law more reflective of the 
profile characteristics exhibited in Figure 3 in order to successfully 
predict the flow field properties occurring in typical turbomachine 
applications . 

Several such modified wall treatments have been proposed 
[8,9,24,25].. Recently Nakayama and Kato [8] derived a wall law, 
including the effects of adverse pressure gradients, from a one- 
dimensional analysis of the turbulent kinetic energy equation with 
gradient diffusion employed to n.^del the near wall shear stress 
variation. This approach was modified by Chen [9] and incorporated as 
an option in the CNS2D computer code. Details of this procedure are 
quite complex and the reader is referred to references [8,9] for a 
complete development. 

In section iv of this report, TAD flow predictions are examined 
using a traditional k-e turbulence model with a conventional wall 
treatment developed by Chen [9]. Comparisons with results of TAD 
analysis using the multiple-scale turbulence model are presented <. id 
phenomenological differences in prediction are displayed. 
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Figure 3 . Turbulent boundary layer profiles for various 
pressure gradients using inner law variables. 
See e.g. reference [28], 
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IV 

RESULTS 


A total of ten CFD analyses of TAD flows within the duct 
configuration depicted in Figure 1 were performed using the CNS2D 
computer code. Three different grid discretizations were employed. A 
coarse 61x21 mesh was initially specified in order to determine 
starting estimates of the flow field characteristics. Near turn detail 
of this mesh configuration is displayed in Figure 5. Longitudinal 
refinement of this discretization was introduced using a 91x21 mesh, 
also displayed in Figure 5. Additional cross-stream detail was 
achieved using a 91x41 mesh. It was found that flow field 
characteristics determined using the refined meshes at significantly 
higher computational overhead differed only slightly (±3%) from 
results obtained using the coarse 61x21 discretization. Therefore, 
flow field properties presented in the balance of this report were 
obtained from computations based on the 61x21 mesh configuration. 

Of the ten CFD analyses performed, seven were conducted to test 
the sensitivity of the turbulence models to variation in model tuning 
parameters (turbulence model constants). Results of these sensitivity 
analyses are presented at the end of this section. The remaining 
three analyses were conducted using specific turbulence model options 
available within the CNS2D computer code. A description of these 
options is given below. 

Option 1 - designated ss k-e, employs a conventional r ingle scale k-e 
turbulence model with the standard logarithmic law used in 
the wall function treatment 

Option 2 - designated mwf, employs a conventional single scale k-e 
turbulence model with a pressure modified wall law used in 
the wall function treatment 

Option 3 - designated ms k-e, employs the multiple scale k-e 
turbulence model described in the previous section together 
with the standard logarithmic law wail function treatment 

The results presented in this report were obtained using the CNS2D 
computer code implementation of these three turbulence model/wall 
function options applied to the TAD configuration displayed in 
Figure 1. 

In all analyses, the working fluid was assumed incomprei : ib_ e 
and the molecular viscosity assumed constant. A fixed Reynolds number 
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value of 10®/ based on the duct width D, was specified. This 
particular value of the duct Reynolds number was selected in order to 
provide computational results for comparison with TAD experimentation 
in progress at Colorado State University (reference NASA contract 
number NAS8-36354) . The target Reynolds number of this experimental 
program is 10 6 . Unfortunately, as of this writing, results at this 
value of the duct Reynolds number were unavailable. 

The experimental water flow facility at Colorado State consists 
of a TAD with a cross section aspect ratio of 10 and with a turn 
centerline radius to duct width ratio, R c /D, equal to 1. The 
apparatus cross section is very similar to the computational 
configuration displayed in Figure 1. A more detailed schematic sketch 
of this test facility is presented in reference {26] . 

A fully developed turbulent duct flow inlet velocity profile, as 
displayed in Figure 4, was assumed for all CFD analyses. A summary of 
the assumed computational boundary conditions is presented below. 

ss k-e and mwf analyses ms k- e analyses 

see Figure 4 
0 _ 3 
k p =k t =1.0 x 10 J 

£p = /2 , e t -e 


The inlet values of k and e in the above table have been normalized by 
U 2 and TP/d respectively. 

Results of the CFD analyses utilizing the various model options 
described above are presented in terms of four dependent variable 
measures and one independent variable phenomenological measure as 
listed below. 

1) Normalized stream function contours (4>* = \1>/J$jb] ) 

2) Normalized pressure contours (P ■ [P-P e ]/lPU z ] # ) 

3) Normalized turbulent kinetic energy contours^tk = k/U^) 

4) Normalized longitudinal velocity profiles (U ■ U/U) 

5) Separation zone normalized size. 

Predicted stream function contours using the various turbulence 
model options are displayed in Figure 6. A most striking 
phenomenological difference is exhibited in the complete absence of a 
downstream recirculation zone for the ss k-e option in Figure 6a. 
Both the multiple-scale model an<] the modified wall function approach 
predict sizeable separation regions near the TAD exit plane. 


Inlet u 
Inlet v 
Inlet k 


see Figure 4 
0 _ 3 
2.0 x 10 J 


Inlet 


C y k 1 * 5 /0.03 


A,. 4-1^4- Un 
wuuxcu « r 


A 
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Separation is observed to occur somewhat over 9 degrees before the 
exit plane is reached using a pressure modified wall law as shown in 
Figure 6b. Using the multiple-scale model, an even longer 
recirculation region is predicted, with flow separation delayed until 
just after the TAD exit. 

The occurrence of a recirculation region, with the potential for 
introducing flow instability and dead zone transport effects, is a 
most important phenomenological event to predict. This dissimilarity 
in the various model predictions of an important physical phenomena 
presents an ideal criterion for evaluation of model capabilities in 
strongly curved internal flows. Application of this criterion must of 
course await reliable experimental data. 

The pressure contours displayed in Figure 7 again highlight 
differences in the model predictions. Pressure recovery downstream of 
the TAD is predicted to occur much faster using the conventional k-e 
model. In all cases, a long pressure recovery region on the order of 
seven duct diameters in length in the exit section is predicted. 
Maximum and minimum pressures, governed almost exclusively by 
centrifugal effects in the turn around section, are virtually 
identical for all the models. 

The turbulent kinetic energy profiles exhibited in Figure 8 
highlight the most striking difference between multiple-scale model 
predictions and results based on the conventional k-e model. This 
difference is not evidenced in the infeed section of the duct. Within 
the turn, however, a substantially reduced value of the turbulent- 
kinetic energy is determined using the multiple-scale model. A much 
reduced turbulent diffusivity is thus indicated when compared with the 
traditional k-e model predictions. Near the duct exit plane, this 
difference is particularly significant, approaching a full order or 
magnitude reduction in v fc . 

Predicted values of the turbulent kinetic energy and diffusivity 
were generally somewhat less in the TAD section with the pressure 
modified wall law than with conventional wall functions as displayed 
in Figures 8a and 8b. The anamolous TAD section k contours are 
strongly affected by sharp mean field velocity gradients, both 
longitudinal and spanwise. Verification of these unusual patterns 
must await additional experimental evidence. 

Mean longitudinal velocity profiles are presented in Figures 9, 
10, and 11 for the TAD entrance, 90 degree line, and exit plane 
respectively. Little difference between model predictions is observed 
at the turn entrance and 90 degree plane, resembling in each case a 
potential flow profile. The multiple-scale model did. however, tend 
to predict lower mean velocities near the concave wall compared to the 
other models. 
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61 x 21 



61 x 21 


Figure 5. Near turn mesh configurations for T^ri analyses 
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** ' Single scale k-e model with standard wall function. 


c) 

Figure 



Single scale k-e model with modified wall function. 



Multiple scale k-e. model with standard wall function. 


7. Normalized pressure contours for TAD analyses, Re=10* 
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a) Single scale k-e model with standard wall function. 




c) Multiple scale k-e model. with standard wall function. 

Figure 8. Normalized turbulent kinetic energy contours for TAD analyses 
Re-106. 
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Considerable variation in the mean field velocity profiles at 
the,exit plane is displayed in Figure 11. The conventional k- model 
profile strongly resembles a fully developed turbulent duct flow 
profile at this point, evidencing an extremely fast adjustment to the 
approached exit section. It is well known that the k equation in the 
standard k-e model exhibits an overly fast response to the mean shear 
stress development [9]. This process is symbiotic, effectively 
accelerating the mean field adjustment as evidenced in Figure 11. The 
profile symmetry exhibited when using the modified wall law with the 
k-e model is also a reflection of this tendency for fast mean field 
adjustment, with exit plane recirculation causing core profile 
overshoot. Only the multiple-scale model results exhibit substantial 
profile assymmetry, indicating a retarded mean field response. 

Figures 12, 13, and 14 display predicted turbulent kinetic 
energy profiles at the turn entrance, 90 degree plane, and exit 
respectively. Inlet profiles were virtually identical except near the 
outer wall where the modified wall function predicted a significant 
profile overshoot when compared with both the single-scale and split 
spectrum models. Previous study [26] of sharply curved TAD flows with 
R c /D = 1 but at reduced flow rates has indicated a tendency of 
computational models to overpredict the turbulent kinetic energy at 
the turn inlet. In this regard multiple-scale model predictions 


cXuibltcu ii O teiiueiiuy 

however, as evidenced 
predicted k values was 
experimental evidence 


foi. downward adjustment at 6=0 degrees, 
in Figure 13, a significant decrease in the 
observed at the mid-turn plane. There is some 
[26] at lower Reynolds numbers to support 


reduced k values near the inner wall. Full k profile reduction as 
predicted by the multiple-scale model has not been observed. Exit 


plane kinetic energy profiles displayed in Figure 14 have the same 
general features as the midplane profiles with the exception of the 


modified wall function profile. The mwf approach provided substantial 
profile overshoot in the exit plane recirculation region. It should 


be observed that all model profiles exhibited in Figures 13 and 14 
differ substantially from previously reported computational results 
and experimental observation [26] at lower Reynolds numbers. 


Predicted pressure profiles at the TAD entrance and 90 degree 
sections were virtually identicax for all the models tested as 
displayed in Figures 15 and 16. These profiles exhibit the 
predominance of a centrifugal ly induced cross stream pressure 
gradient. At the exit plane, use of the modified wall function model 
significantly retarded the inner wall pressure recovery while 
enhancing outer wall pressure relaxation as shown in Figure 17. 
Similar remarks apply to results based on the multiple-scale model as 
compared to the conventional k-e model predictions exhibited in Figure 
17. Further experimental evidence at high Reynold® numbers will l >e 
required in order to ascertain the validity of any of hese 
significantly differing cross stream profiles. 
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The computational results described above highlight two striking 
differences in the various model results. The first is the downstream 
separation-reattachment region recovered using both the modified wall 
function approach and the split spectrum model which is absent from 
results based on a conventional k-e model treatment. Although there 

. , • ±_ - 4 4-a Ar vaf nf 9 th9 SXlSt^HC? O I SUCH 

is lime tuitcia ww — rr-" - - V /• ** 

a zone at high flow rates, the computational differences do present an 
effective criterion for comparative model evaluation when experimental 
data becomes available. 

The second major difference involves the substantial reduction 
in turbulent diffusivity associated with use of the multiple-scale 
model. There exists some evidence obtained at lower Reynolds numbers 
to support a reduction in the turbulent kinetic energy and diffusivity 
in sharply curved flows [26]. Additional information is required in 
order to make the appropriate quantitative comparisons needed for 
model verification. 

In addition to the basic model comparisons described above, 
several computational analyses were performed in order to test the 
sensitivity of the multiple-scale approach to the model tuning 

parameters, o k , o e , Cp 1# C p2 , C t l' and c t 2 ‘ 4 ^ h ® ][* 
these parameters in turn was increased 3\ ±n ordei. cw 
effects on computational predictions. Except for C pl , little 
variation either in the dependent variable measures or in the 
phenomenological recirculation zone measure was observed. A 261 
increase in the length of the recirculation zone was obser «pon 
changing C n1 3% as shown in Figure 18. Values of the dependent 
variable wire little modified anywhere except near the separation 
region. This is indicative of an inherent instability associateu wiv*. 
this region. The relative insensitivity °. f ,^he muitipre-scaie 
approach to other parameter variation is an indication oj. a iOuusi. 
model. Little objective evidence is available to confirm accuracy for 
the flow conditions tested. As with most turbulent flow computation, 
more experimental verification is necessary. 
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Figure 9. TAD inlet velocity profiles, Re=10 6 . 



Figure 10. TAD 90° velocity profiles, Re=10 6 . 
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OBSERVATIONS AND RECOMMENDATIONS 


1) The CNS2D computer code provides a flexible 
computational tool for BFC coordinate analysis. 


and robust 


2) At a Reynolds number of 10° both the multiple-scale k-t model and 


6 

VV W1A WAAW — - § 

the conventional k-e model with a pressure modif led wall 
treatment provide TAD predictions which differ significantly from a 
conventional k-e model approach using a logarithmic law wall 
Both modifications of the standard two equation model predict flow 
seoaration at the TAD exit plane, even at this elevated Re value. 
This is in contrast to the no separation prediction of the standard 

model . 

3) The multiple-scale turbulence model predicts significantly reduced 
turbulent diffusivity in sharply curved flow regions at high Re 

values . 

4) Both the multiple-scale k-e model and the conventional k-e model 
with pressure modified wall function treatment provide viable 
turbulence modeling modifications and the potential for ®^ lflcant 
model improvement with little additional computational overhead. 


2GG i t 1 ‘Jlld i i-- - 

should be tested. 


5) Modification or me 

terms due to extra, curvature induced, strain 

6) A modified wall function approach appropriate to the multiple-scale 
k-e model should be developed in order to simultaneously see, ount for 
spectral non-equilibrium and adverse pressure gradient in sharply 

curved flows. 

7) Additional TAD experimentation is needed in order t0 pr 

basis for testing the applicability of improved turbulence models to 
sharply curved flows at very high Reynolds numbers. 
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